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I.  SUMMARY 


This  is  a  study  of  some  fundamental  aspects  of  the  structural 
dynamics  and  vibration  control  of  large  structures.  One  focus  is  the 
development  of  a  limiting  performance  formulation  with  minimum  settling 
time  which  can  accept  multiple  design  objectives  efficiently.  This  new 
formulation  is  intended  to  meet  the  need  of  rather  comprehensive  design 
objectives  for  the  control  of  large  space  structures.  Another  objective 
of  the  study  is  to  develop  a  systematic  way  of  designing  a  control 
system  based  on  the  limiting  performance  characteristics.  An  indirect 
synthesis  methodology  is  proposed.  It  is  shown  that  closed  loop  control 
laws  can  be  based  on  the  optimal  response  trajectories  in  the  time 
domain.  The  method  is  successfully  applied  to  the  control  of  proof-mass 


actuators. 


II.  RESEARCH  OBJECTIVES 


1 .  A  Comprehensive  Limiting  Performance  Study 

A  comprehensive  limiting  performance  study  was  to  be  performed.  In 
order  for  the  varied  and  stringent  design  objectives  of  large  space 
structure  control  problems  to  be  met  in  terms  of  limiting  performance, 
it  was  considered  essential  that  a  broader  limiting  performance 
formulation  be  developed.  In  addition  to  the  capability  of  taking 
various  objectives  the  resulting  trajectories  were  required  to  be 
unique.  In  conjunction  with  the  synthesis  study  of  feedback  control 
laws,  limiting  performance  characteristics  in  the  modal  coordinates  were 
to  be  provided. 

2.  Synthesis  of  Feedback  Control  Laws 

Since  the  limiting  performance  characteristics  with  minimum 
settling  time  could  give  optimal  response  trajectories  in  the  time 
domain  even  when  the  problems  included  constraints,  it  was  anticipated 
that  a  methodology  to  give  control  laws  based  on  these  characteristics 
would  be  an  efficient  design  tool.  An  objective  of  this  research  was  to 
develop  such  a  methodology.  In  order  to  assure  robustness  and 
practicality,  it  was  felt  that  such  control  laws  should  be  of  closed 
loop  form. 


2 


3. 


The  control  laws  developed  were  to  be  applied  to  the  control  of 


proof-mass  actuators.  Since  the  control  characteristics  of  proof-mass 


actuators  are  dominated  by  the  constraints  on  the  control  force  and  the 


rattlespace,  the  control  laws  should  be  able  to  be  designed  to  handle 


the  constraints  efficiently. 


The  limiting  performance  is  defined  as  the  minimum  peak  value  of 
certain  responses  while  other  system  responses  are  constrained.  Min-max 
norms  are  utilized  in  the  performance  index  to  achieve  the  limiting 
performance  characteristics.  Since  the  min-max  norm  employed  to 
minimize  the  maximum  peak  value  among  the  peak  values  of  the  selected 
responses  does  not  provide  unique  peak  values  except  the  maximum  peak 
value,  an  norm  with  weights  is  proposed  here  to  avoid  the 
non-uniqueness  of  the  peak  values  and  to  provide  flexibility  in 
selecting  the  performance  index.  Also,  it  is  shown  here  how  to 
supplement  min-max  norms  with  additional  performance  measures  to  settle 
the  responses  in  minimum  time.  The  formulations  presented  provide  a 
general  approach  to  find  unique  optimal  response  trajectories.  The 
resulting  characteristics  are  referred  to  as  limiting-performance/ 
minimum-settling-time  (LP/MST)  characteristics. 


Problem  Statement 


A  linear  vibrating  system  can  be  represented  by  a  system  of  second 


order  linear  differential  equations 


where  x  is  an  n-dimensional  vector  and  M,  C,  and  K  are  mass,  damping, 
and  stiffness  matrices  with  appropriate  dimensions,  f.(t)  is  tin 
nf-dimensional  vector  representing  a  forcing  function,  and  F  is  a 
coefficient  matrix  which  is  used  to  place  the  forcing  function  in  the 
equations  of  motion.  By  replacing  portions  of  the  system  to  be  designed 
by  generic  or  control  forces  which  can  represent  tiny  configuration,  the 
limiting  performance  characteristics  of  the  system  can  be  found.  With 
the  control  forces,  u(t).  the  system  to  be  studied  can  be  described  by 
the  equations  of  motion 


Mx  +  Qc  +  Kx  +  Vu(t)  =  Ff ( t)  (2) 

where  u(t)  is  an  nu-dimensional  control  force  vector,  V  is  a  coefficient 
matrix,  and  C  and  K  represent  the  resulting  damping  and  stiffness 
matrices  after  replacing  portions  of  Eq.  (1)  with  control  forces. 

Define  a  state  vector 


T  -T  T 
s  =  [x  X  J 


Then,  Eq.  (2)  can  be  represented  by  the  state  equations 


s(t)  =  As(t)  +  Bu(t)  +  Cf(t) 


(4) 


where  s(t)  is  a  2n-dimensional  state  vector.  A,  B,  and  C  are  constant 
coefficient  matrices  represented  by 


where  I  is  the  identity  matrix  and  0  is  the  null  matrix. 

The  problem  is  to  find  an  optimal  control  u  (t)  based  on  a 
performance  index  as  the  weighted  sum  of  p  system  responses 


minimize 


J  =  wj  =  WjJj  +  w2J, 


+  w  J 
P  P 


where 


1  =  [  Jj  J2  . • •  Jp  ]T  =  tQ  <  t  <  tf  IPjS  +  PgU  +  P3f |  (9) 

w  =  [  Wj  w2  ...  wp  ]T  (10) 

where  t  and  t^.  are  given  initial  and  final  times,  P^ ,  P2  and  P3  are 
prescribed  constant  coefficient  matrices,  and  w  is  a  prescribed  constant 
coefficient  vector  containing  weights.  Equations  (8)  and  (9)  represent 
and  min-max  norms,  respectively.  The  solution  is  to  be  computed  such 
that  the  peak  values  of  j  are  unique  and  the  response  trajectories  after 
the  peak  values  of  responses  are  settled  in  minimum  time. 


Constraints  are  imposed  on  the  dynamic  system  under  study.  The 


format  of  the  constraints  is 


A 


( 


XL  <  QjS  +  QgU  +  Q3£  <3^  for  tQ  <  t  <  tf  (11) 

where  y^  and  y^  are  nc-dimensional  lower  and  upper  constraint  vectors 
and  ,  Qg  and  are  constant  coefficient  matrices.  The  quantity  nc  is 
the  number  of  constraints 

Linear  Programming  Formulation 

Since  the  problem  defined  above  has  a  linear  performance  index  and 
linear  constraints,  standard  linear  programming  can  be  used  to  solve  the 
optimization  problem,  provided  that  the  problem  has  not  been  stated  to 
be  physically  over-constrained.  In  this  section,  a  linear  programming 
formulation  without  imposing  minimum  settling  time  is  studied.  If  the 
system  in  Eq.  (4)  is  discretized  using  N  uniform  time  intervals,  a  set 
of  state  difference  equations  is  obtained  for  k  =  1,2 . N 

s(k)  =  Gs(k-l)  +  H[Bu(k)  +  Cf(k)]  (12) 


where 


s(k)  =  state  vector  at  time  t^  =  kh 

u(k) ,  £(k)  =  control  and  external  excitation  vectors,  assumed  to  be 


constant  over  the  interval  t^_j  <  t  <  t^ 


Ah 


G  =  e 
H  =  /h  eA(h_T)dT 

J  A 


h  =  time  interval  =  t^  -  t^_j 


"wxwuv'JTfvr.Truv'.v'.v. 


iv 


The  state  vector,  at  any  time  t  =  t^.  can  be  expressed  as  a  function  of 

the  initial  state  s(0),  the  control  history  u(l),  u(2) . u(N).  and 

the  external  excitation  £(1),  £(2),  ....  £(N).  For  k  =  1,2 . N 

s(k)  =  Gk  s(0)  +  2  Gk-j  H[Bu( j)  +  C£(j)]  (13) 

j=l 

The  constraints  in  Eq.  (11)  are  discretized  similarly 

XL(k)  <  QjKk)  +  Qgufk)  +  Q3£(k)  <  y^k)  for  k  =  1.2 . N  (14) 

Having  discretized  the  state  equations,  these  constraints  can  only  be 
enforced  at  the  end  of  each  interval.  At  "intermediate"  times,  they  may 
be  violated. 

Consider  now  the  linear  programming  formulation  of  the  limiting 
performance  problem.  The  objective  function  of  Eq.  (9),  which  reflects 
the  min-max  norm,  can  be  converted  into  a  constraint  set.  Since  j  is 
the  vector  of  maximum  values  of  |PjS(k)  +  Pgufk)  +  P^£(k) |  for  all  k, 

|PjS(k)  +  P^fk)  +  P3£(k)|  <  j  (15) 

To  place  the  norm  of  the  performance  index  in  Eq.  (8)  into  linear 
programming  form,  define  a  (p+nc)  x  p  augmented  matrix  J  which  has 
elements  of  vector  j  along  its  diagonal  as  the  only  non-zero  elements. 


V 


J  = 


• 

.  0 

1 

1 

o  I 

1 

1 

1 

1  o 

1 

1 

<  <- 

0  . 

o 

■  o 

•  o 

(16) 


Also,  to  merge  constraints  in  Eqs.  (14)  and  (15),  define  augmented 
matrices 


Q1  = 


[«;]• 


(17) 


and  augmented  vectors 


^(k)  =  [QT  yJ(k)]T  and  iL(k) 


[0T  xJ(k)]T 


(18) 


where  0  is  a  p-dimensional  null  vector.  Then  Eqs.  (14)  and  (15)  can  be 
merged  as 


-J  +  Qxs(k)  +  QgUfk)  +  Q3i(k)  <  %(k) 
-j  -  Q^k)  -  Qgufk)  -  Q3f(k)  <  -iL(k) 


(19) 

(20) 


Substituting  Eq.  (13)  into  Eqs.  (19)  and  (20),  the  constraint  equations 
have  the  control  sequence  u(k)  as  the  only  unknowns.  For  k  =  1,2 . N 


~k  ji 


-J  +  Q,  2  G  JHBu( j)  +  CLu(k)  < 

j=l 

,  k  ,  . 

iyCk)  -  QjG^sfO)  -  Qj  2  GK_JHC£(j)  -  Q3£(k) 

-  -  k  k-j  - 

-J  -  Qj  2  Gk  JHBu(j)  -  Qgufk)  < 

-iL(k)  +  Q^sfO)  +  Qj  ^Gk_JHC£(j)  +  Q3£(k) 


(21) 


(22) 


Since  the  control  u(k)  should  be  unrestricted  in  sign,  a  standard  linear 
programming  technique  for  dealing  with  such  variables  is  introduced. 

Def ine 


u(k)  =  U  (k)  -  U  (k) 


(23) 


where 


U  (k)  >0  and  U" (k)  >  0 


(24) 


To  place  this  optimization  problem  into  a  standard  linear  programming 
form,  define 


-  ii) 


(25) 


where 


u  =  [U+(1)T  u  (i)T  u+(2)T  u  (2)T  ...  u+(n)t  u  (n)t]t  (26) 


and 


T  r  T  nT  T 
c  =  [  w  0  ] 


(27) 
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where  0  is  a  2N-dimensional  null  vector.  The  vector  w  can  be  selected 

to  provide  unique  values  of  peaks,  i.e.,  J^,  Jg . Jp -  Then  the  linear 

programming  problem  is  to  minimize 

J  =  cTz  (28) 

subject  to  the  constraints 

Hz  <  b  (29) 

where  H  and  b  are  a  2N*(p+nc)  x  (2N*nu+p)  coefficient  matrix  and  a 
2N*(p+nc)  coefficient  vector,  respectively,  representing  constraints  of 
Eqs.  (21)  and  (22). 

Limiting-Performance/Minimum-Settling-Time  Formulation 

The  unique  peak  values  can  be  obtained  based  on  the  linear 
programming  formulation  of  the  previous  section.  Formulations  are 
studied  to  achieve  the  minimum  settling  time  characteristics  after  the 
peak  values.  Two  methods  are  available  to  achieve  the  minimum  settling 
time  and  they  can  produce  identical  minimum  settling  time  responses. 

One  utilizes  additional  performance  indices  and  the  other  utilizes 
additional  constraints.  The  additional  performance  indices  and 
constraints  can  be  employed  independently  or  simultaneously  depending  on 
the  requirements  of  a  problem.  However,  in  order  to  minimize  the  size 
of  the  linear  programming  problem,  additional  performance  indices  for 
the  responses  selected  in  the  performance  index  and  additional 


?TT 


ilA] 
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constraints  for  the  constrained  variables  are  suggested  to  be  employed. 

To  provide  the  minimum  settling  time  characteristics  for  the 
responses  in  the  performance  index,  the  performance  index  of  Eq.  (8)  is 
separated  into  two  performance  indices.  Assume  that  the  peak  values  of 
the  limiting  performance  occur  during  the  time  between  t  and  tt.  One 
of  the  performance  indices,  referred  to  as  the  transient  performance 
index,  is  given  by 


rt  T.t 
J  =  w  j 


where 


^  =  t  <mJX<  t  |PjS  +  PgU  +  P3£ | 


where  t^  is  the  time  limit  for  the  transient  period.  The  other, 
referred  to  as  the  steady-state  performance  index,  is  defined  as 


Ts  T  .s 
J  =  w  j 


where 


=  t  <"^X<  t  | PjS  +  PgU  +  P3f | 


Note  that  the  same  w  and  P^ ,  Pg,  P3  are  used  in  the  performance 
indices.  The  transient  performance  index  is  used  to  give  unique 
response  trajectories  up  to  the  point  of  the  peak  values  and  the 
steady-state  performance  index  is  used  to  settle  the  response 
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trajectories  after  the  peaks  within  the  desired  values  in  minimum  time. 
See  Fig.  1. 


Fig.  1  A  Response  Trajectory  with  Transient  and  Steady-State 
Performance  Indices  (i  =  l,2,...,p) 


Now.  the  global  performance  index  is  defined  by 


jg  =  j‘  ♦  js 


The  objective  functions  of  Eqs.  (31)  and  (33),  which  reflect  min-max 
norms,  are  also  entered  into  the  constraints.  Since  is  the  vector  of 
the  maximum  values  of  |PjS  +  PgU  +  P^f |  for  tQ  <  t  <  tt  and  so  is  jS  for 


IV I 

I 


t  t  <  tj.,  the  constraints  from  the  objective  functions  are,  if  they 
are  discretized  using  the  same  discretization  of  Eq.  (12), 

|PjS(k)  +  P^k)  +  P3f(k)|  <  for  k  =  0.1 . Nt_1 

and  (35) 

IP^k)  +  P^k)  +  P3f(k)|  <  js  for  k  =  Nt,Nt+l . N 

where  NC  is  the  discretized  time  limit  for  the  transient  period. 


Fig.  2  A  Response  Trajectory  with  Transient 
and  Steady-State  Constraints 


To  achieve  minimum  settling  time  characteristics  for  the  responses 
in  the  constraints,  consider  two  sets  of  constraints.  One  set  of  them, 
referred  to  as  the  transient  constraints,  is  used  to  represent  the 
original  constraints  in  Eq.  (11)  from  tQ  to  t^.  An  additional 
constraint  set,  referred  to  as  steady-state  constraints,  is  imposed  from 
t^  to  t. 


*SL  1  Ql“  +  V  +  Q3L  1  *SU  for  \  1  t  1  *f 


(36) 


where  Ygj^.  are  nc-dimensional  coefficient  vectors  representing  lower 
and  upper  bounds  of  steady-state  constraints.  The  steady-state 
constraints  represent  the  desired  values  of  the  responses  after  the 
peaks.  The  responses  are  required  to  be  settled  within  these  values  in 
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the  smallest  value  of  t 


Note  that  the  matrices  , 

(36)  are  identical  with  those  in  Eq.  (11).  See  Fig.  2.  for  an 
illustration. 


Qg.  and  Q3  in  Eq. 


To  take  the  new  performance  index  of  Eq.  (34)  into  linear 
programming,  define  p  x  p  matrices  J  and  J  which  have  elements  of 
vectors  ^  and  ^  ,  respectively,  along  their  diagonals  as  only  non-zero 

A  £  Ag 

elements.  Also,  define  (p+nc)  x  2p  matrices  J  and  J  such  as 


J*  =  [  jt  °  1  (37) 

L  0  0-1 

5s  =  [  0  JS1  (38) 

L  0  0-1 


where  0  represent  null  matrices.  The  two  sets  of  constraints  are 
discretized  and  merged  into  the  constaints  from  the  performance  index, 
Eq.  (35).  using  augmented  matrices  and  vectors  of  Eqs.  (17)  and  (18). 
Then,  for  k  =  1,2,...,^  the  constraints  are 


-?  +  Q^k)  +  QgUfk)  +  Q3f (k)  <  ^(k)  (39) 

-J*  -  Qxs(k)  -  Q^k)  -  Q3f(k)  <  -yL(k)  (40) 

and,  for  k  =  Nt.Nt+l . N,  the  constraints  are 

-JS  +  QjS(k)  +  Q^k)  +  Q3f(k)  <  ^j(k)  (41) 

-JS  "  Q^fk)  -  QguCk)  -  Q3f  (k)  <  -igL(k)  (42) 
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Substituting  Eq.  (13)  into  Eq.  (39)  to  Eq.  (42),  the  constraint 
equations  have  control  sequence  u(k)  as  only  unknowns. 


To  place  this  optimization  problem  into  a  standard  linear 
programming  form,  define 


z  = 


(43) 


where  u  is  given  by  Eq.  (26)  and 


~r 

c 


r  T  T  -.T-, 
=  [  w  w  0  ] 


(44) 


where  0  is  a  2N-dimensional  null  vector.  Then  the  linear  programming 
problem  becomes:  minimize 


(45) 


subject  to  the  constraints 


Hz  <  b  (46) 


where  H  and  b  are  a  2N*(p+nc)  x  (2N*nu+2p)  coefficient  matrix  and  a 
2N*(p+nc)  coefficient  vector,  respectively,  representing  constraints  of 
Eq.  (39)  to  Eq.  (42). 
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Limiting  Performance  Based  Feedback  Control  Synthesis  and  Its 


Application  to  the  Control  of  Proof-Mass  Actuators 


The  control  characteristics  of  proof-mass  actuators  are  dominated 
by  the  constraints  on  the  control  force  and  the  rattlespace.  This 
constrained  control  problem  raises  a  serious  difficulty  in  designing  a 
system  using  conventional  control  system  design  methods  such  as  linear 
control  or  optimal  control.  It  has  been  known  that  the  limiting 
performance  formulation  with  the  minimum  settling  time  can  provide  an 
optimal  open  loop  control  law  for  the  proof-mass  actuators.  In  order  to 
enhance  the  practicality  and  robustness  of  the  open  loop  control  and  to 
provide  a  systematic  approach  to  handle  constrained  control  problems, 
feedback  control  synthesis  methods  based  on  the  limiting  performance 
trajectories  are  developed.  Some  control  laws  developed  with  these  new 
approaches  for  the  control  of  proof-mass  actuators  are  presented. 


Suboptimal  Design  Synthesis  of  Feedback  Control  Systems 


Methods  of  designing  suboptimal  feedback  control  laws  based  on  the 
LP/MST  characteristics  are  presented  here.  Assume  that  the  system  is  to 
be  controlled  by  active  actuators  and  the  actuators  are  subject  to 
constraints.  The  active  controllers  are  assumed  to  be  controlled  by  a 
constant  state  variable  feedback  and  the  feedback  gain  matrix  is  found 
based  on  the  LP/MST  characteristics.  Two  distinct  design  synthesis 
methods  to  construct  feedback  control  laws  are  studied.  For  general 
equations  of  motion  including  systems  with  a  non-proportional  damping 
matrix,  a  suboptimal  feedback  gain  matrix  is  found  using  the  constrained 


curve  fitting  technique  based  on  approximating  optimal  response 
trajectories.  For  decoupled  equations  of  motion,  e.g.,  systems  with  a 
proportional  damping  matrix,  another  synthesis  method  can  also  be 
utilized.  In  this  method,  closed  loop  eigenvalues  are  first  extracted 
using  the  curve  fitting  based  on  approximating  optimal  response 
trajectories  and  the  eigenvalues  are  assigned  using  an  eigenvalue 
assignment  method  based  on  a  constant  state  variable  feedback  to  find 
the  feedback  gain  matrix. 

Identification  of  Feedback  Control  Based  on  the  Optimal  Response 
Trajectories 

One  of  the  two  methods  to  construct  feedback  control  laws  based  on 
the  LP/MST  characteristics  is  studied  here.  With  this  method,  the 
constant  state  variable  feedback  gains  are  obtained  directly  by  fitting 
the  optimal  response  trajectories.  This  method  can  handle  equations  of 
motion  with  a  non-proportional  damping  matrix.  However,  there  is  no 
guarantee  that  the  resulting  system  will  be  stable. 

Computation  of  LP/MST  Characteristics 

The  first  step  is  to  compute  the  LP/MST  characteristics.  By 
replacing  portions  of  a  system  to  be  designed  by  generic  forces  which 
can  represent  any  configuration,  the  LP/MST  characteristics  of  the 
system  can  be  found  as  described  previously. 


Identification  of  Constant  Feedback  Gain  Matrix 

The  optimal  LP/MST  trajectories  s*(k)  =  [x*(k)^  x*(k)^]^  and  u*(k) 


(k  =  1,2 . N)  are  used  as  a  starting  point.  Assume  that  the  active 

controllers  are  controlled  by  a  constant  state  variable  feedback,  i.e., 

u(k)  =  Ds(k)  (47) 

where  D  is  an  nu  x  2n  constant  feedback  gain  matrix.  Since  nu 
controllers  are  considered  and  optimal  control  forces  are  available  for 
each  controller,  it  is  possible  to  proceed  controller  by  controller. 

For  controller  j  (j  =  l,2,...,nu),  a  suboptimal  linear  control  law  to  be 
determined  is  described  by 

Umax.  for  u®(k)  >  Umax. 

J  Jv  J  j 

Uj(k)  for  Uminj  <  u®(k)  1  U,naxj  (48) 

Uminj  for  u^(k)  <  U|minj 

where 

«j(k)  =  djS^k)  (49) 

dj  =  row  of  D  matrix 

Note  that  Eq.  (48)  serves  to  enforce  a  constraint  on  the  control  force 
of  controller  j 

Uminj  <  Uj(k)  <  UmaXj  (50) 

Then  the  problem  becomes  to  find  the  optimal  feedback  gains  d.  (for  j  = 

J 

1,2 . nu)  which  minimize 


where 


k 


P 


Aj(k)  =  |u*(k)  -  u*(k) | 


(52) 


To  find  dj  efficiently  a  curve  fitting  method  based  on  a  least  squares 


residual  function  of  Eq.  (51)  can  be  used. 


Identification  of  Feedback  Control  Based  on  the  Assignment  of  the 
Optimal  Closed  Loon  Eigenvalues 


In  this  section  it  is  assumed  that  the  equations  of  motion  can  be 
decoupled,  i.e.,  the  damping  matrix  is  proportional  or  the  system  is 
undamped.  In  this  method,  the  optimal  response  trajectories  are 
calculated  in  the  modal  coordinates  and  the  curve  fitting  technique  is 
used  to  find  closed  loop  eigenvalues  for  each  mode.  Then,  the  state 
variable  feedback  gains  are  obtained  using  an  eigenvalue  assignment 
technique.  This  method  can  guarantee  the  stability  of  the  resulting 
gains.  Also,  the  constraints  on  the  response  variables  can  be  enforced 
during  the  design  process. 


Computation  of  LP/MST  Characteristics  in  Modal  Coordinates 
Replace  portions  of  a  system  to  be  designed  by  generic  forces, 
decouple  the  equations  of  motion  by  using  the  normal  mode  method,  and 
compute  the  limiting  performance  characteristics  in  the  modal 
coordinates . 
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where 


< 


A  = 

A 

B  = 

A 

c  = 


0  I 

~_1  ~ 

-M  lK  -M  C 

-H-‘V  ] 

-i-'f  ] 


I 


(59)  ; 

I 

l 

1 

t 

(60) 

I 

i 

t 

(6i)  : 


The  performance  index  and  the  constraints  can  also  be  transformed  with 
the  relationship 


A  A 

s(t)  =  <J>s(t) 


(62) 


I 


Then  the  problem  can  be  solved  in  the  modal  coordinates  using  one  of  the 
formulations  given  previously. 


►  Identification  of  Suboptimal  Closed  Loop  Eigenvalues 

The  optimal  limiting  performance  trajectory  for  each  mode  in  modal 

coordinates.  qr(k)  (r  =  1,2 . n;  k  =  1,2 . N),  is  used  as  a  starting 

point.  The  first  task  is  to  find  optimal  f  and  w  for  each  mode  which 


will  minimize 


where  the  measure  of  approximation  is 


A;(k)  =  |q*(k)  -  q;(k) I 


where 


q;(k)  =  q;{cr  cr} 

=  The  analytical  solutions  for  the  time  response  of  the 
decoupled  SDOF  system. 


To  find  optimal  f  and  a  least  square  curve  fitting  technique  is 
employed.  Now  the  closed  loop  eigenvalues  for  each  mode  can  be 
calculated  by  the  characteristic  equation 


*2  „  *  w  w  *2  ,  _ 

X  +  2f  u  X  +  <j  =0,  r  =  1,2, 
r  r  r  r  r 


_*  *  *  /,  _* 
=  -{,<i)  +  l  u  v  1  -  f 


Since  the  closed  loop  eigenvalues  are  now  known  for  each  mode, 
eigenvalue  assignment  techniques  can  be  used  to  find  the  feedback  gain 
matrix. 


Assignment  of  Closed  Loon  Eigenvalues 

Several  assignment  techniques  are  available,  often  differing 
according  to  the  characteristics  of  eigenvalues  to  be  assigned.  From 
Eq.  (4)  the  system  to  be  designed  is  written  as 


s(t)  =  As(t)  +  Bu(t)  +  Cf(t) 


(67) 


and  the  problem  is  to  find  the  feedback  gain  matrix  D  such  that 

u(t)  =  Ds(t)  (68) 

which  will  produce  the  specified  closed  loop  eigenvalues  (r  =  1,2,... 
,n)  of  Eq.  (66)  from  the  matrix 


As  a  proof-of-the-concept  example,  the  suboptimal  design  synthesis 
method  is  employed  to  find  suboptimal  feedback  control  laws  for  the 
proof-mass  actuators. 


Proof-Mass  Actuators 


One  of  several  designs  for  a  proof-mass  actuator  is  shown  in  Fig. 
3.  It  incorporates  two  doughnut  shaped  samarium  cobalt  magnets  and  an 
annular  soft  iron  yoke  which  combines  the  functions  of  moving  mass  and 
magnet.  This  assembly  moves  over  a  fixed  coil  whose  current  is 
controlled  to  produce  the  required  force  on  the  magnet,  and  thus  the 
required  reaction  on  the  structure  to  be  damped.  Given  this  coil  size, 
the  maximum  force  available  is  limited  by  the  current  which  can  be 
carried  by  the  coil  without  overheating.  In  the  prototype  design,  it 
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Problem  Statement 

The  system  to  be  considered,  a  single  DOF  mass-spring  system,  is 
shown  in  Fig.  4.  A  proof -mass  actuator  is  attached  at  the  end  of  the 
mass.  The  equations  of  motion  of  the  system  in  Fig.  4  are 

tf 

Mxj  +  KXj  =  u  (70) 

** 

m^  =  -u 

The  natural  frequency  of  the  spring-mass  system  is  chosen  to  be 


w  =  V  K/M  =  2H  rad/s  =  1  Hz 


(71) 


M  =  2.78  kg 
m  =  0.278  kg 


Fig.  4  A  Single  DOF  Mass-spring  System  with  the  Proof-mass  Actuator 


Let  d  =  Xg  -  x^  and  introduce  the  state  vector  s  -  [x^  Xj  d  d]  . 
set  of  first  order  differential  equations  is  obtained. 


s  =  As  +  Bu 


(72) 
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where 


A  = 


'  0 

1 

0 

0  ‘ 

0 

2 

0 

0 

0 

B  = 

1/M 

0 

0 

0 

1 

0 

2 

(J 

0 

0 

0 

M+m 

L  Mm  J 

(73) 


Due  to  the  physical  limitations,  the  proof-mass  actuator  has  constraints 
on  the  distance  it  can  travel  and  the  force  that  can  be  generated. 


I  d  I  <  d  =  1.27  x  10  2m 
1  1  -  max 

I  u  I  <  u  =  2  N 
1  1  ~  max 


The  initial  conditions  are  given  by 


s(0)  =  [0.01m  0  0  0]J 


(74) 


(75) 


The  problem  is  described  in  two  phases.  The  first  phase  is  to  find 
the  optimal  trajectories  of  control,  u  (t),  and  the  state  variables, 
s  (t).  to  bring  the  displacement  of  the  mass  M,  x^,  within  2%  of  the 
initial  value  in  the  minimum  settling  time  and  the  second  phase  is  to 
find  suboptimal  feedback  control  laws  based  on  the  optimal  trajectories 
using  the  methods  for  suboptimal  design  synthesis  of  feedback  control 
systems . 

Computation  of  LP/MST  Trajectories 

The  computation  can  be  done  using  one  of  the  formulations  given 
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on  fitting  optimal  response  trajectories  is  used.  Since  nu  =  1  and  2n 
4,  the  problem  becomes  to  find  a  vector  d  such  as 


u(k)  =  d  s(k) 


(76) 
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Using  a  curve  fitting  technique,  the  suboptimal  feedback  gains  are 
obtained  as 

d  =  [  -81.86  11.50  30.26  8.017  ]  (77) 

and  the  resuting  closed  loop  eigenvalues  are 

-107.4,  -4.821.  -0.5277+5. 956 i  (78) 

The  constraint  on  the  control  force  is  enforced  by  Eq.  (48)  with  Eq. 
(74).  The  time  response  using  these  feedback  gains  for  given  initial 
conditions.  Eq.  (75),  is  shown  in  Fig.  6.  Note  that  the  constraint  on 
the  relative  displacement  is  violated. 


0.0  Time  [sec]  6.0 

Fig.  6  Time  Responses  Using  Suboptimal  Feedback  Gains 
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To  consider  the  influence  of  knocking  between  the  proof-mass  and  the 
housing,  time  response  is  simulated  by  putting  elastic  impacts  when  the 
proof-mass  violates  the  limit  on  the  rattlespace  and  is  shown  in  Figs.  7 
and  8.  Note  that  the  knocking  does  not  cause  instability. 


Identification  of  Feedback  Control  Based  on  the  Assignment  of  Optimal 
Closed  Loop  Eigenvalues 

The  LP/MST  trajectories  shown  in  Fig.  5  are  employed  to  find 
suboptimal  feedback  gain  matrix  D.  In  this  subsection  the  method  based 
on  the  assignment  of  optimal  closed  loop  eigenvalues  is  used.  Since  the 
equations  of  motion  are  decoupled,  the  normal  mode  method  is  unnecessary 
for  this  problem.  Optimal  damping  ratios  and  natural  frequencies  for 
the  mass-spring  system  and  the  proof-mass  are  obtained  based  on  the 

&  K  X 

measure  of  approximation,  Eq.  (64).  The  resulting  Q  and  u  are  = 
0.1087  and  =  6.3641  for  the  mass-spring  system  and  fg  =  0.0381  and 
=  6.1479  for  the  proof-mass  actuator.  Now,  the  closed  loop  eigenvalues 
to  be  assigned  are  determined  from  Eq.  (66) 


> 


9 


\  2  =  ~°-6918  ±  6.326i 


A-  .  =  -0.2342  +  6. 143i 
3.4 


(79) 


Using  an  eigenvalue  assignment  technique,  e.g.,  a  method  based  on  the 
generalized  control  canonical  form,  a  feedback  gain  matrix  D  of  Eq.  (68) 
can  be  obtained  as 
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D  =  [  9.827  0.3844  11.605  0.5517  ] 


(80 


and  the  resulting  time  responses  are  shown  in  Fig.  9.  Note  that  the 
constraint  on  the  rattlespace  is  not  violated.  However,  more  time  is 
required  to  damp  out  the  disturbance  compared  to  the  previous  time 
responses  of  Fig.  7. 


Control.  u(t) 


Displacement.  Xj(t) 

WM \AMam 


II  V  V  '■> 


^M/VVWvwwwn, 


Rattlespace.  d(t) 


Time  [sec] 


Fig.  9  Time  Responses  Using  Suboptimal  Feedback  Gains 


Position  Lood  Based  Desien  Synthesis  of  Feedback  Control 


The  two  methods  shown  previously  give  systematic  approaches  for 
designing  feedback  control  systems  based  on  the  LP/MST  trajectories. 
Although  the  gains  identified  by  fitting  the  optimal  response 
trajectories  give  a  better  damped  result,  the  constraints  were  violated 


V 


7. 


77P7F, 


and  stability  would  not  be  guaranteed  for  any  control  problems  using 
this  method.  On  the  other  hand,  the  gains  obtained  by  assigning  optimal 
closed  loop  eigenvalues  give  a  stable  control  law  without  violating 
constraints,  but  more  time  is  required  to  damp  out  the  disturbance.  To 
achieve  stable,  well-damped  performance  without  violating  constraints, 
the  idea  of  a  position  loop  is  considered  in  conjunction  with  the 
identification  method  based  on  fitting  optimal  response  trajectories. 


Mass-Spring  System 
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Fig.  10  Block  Diagram  of  the  System 
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Fig.  11  Block  Diagram  with  a  Position  Loop 


The  equations  of  motion  in  Eq.  (70)  can  be  represented  by  a  block 
diagram.  As  implied  in  Fig.  5,  the  optimal  motion  of  the  proof-mass 
along  its  finite  length  track  is  rather  simple  compared  to  the 
trajectory  of  the  control,  u(t).  Therefore,  instead  of  utilizing  the 
control  force  as  an  input  to  the  system,  the  relative  displacement  of 
the  proof-mass  is  used  as  the  input  control  variable.  The  block  diagram 
of  Fig.  10  is  expanded  by  attaching  a  position  loop  to  the  plant,  as 
described  by  the  dashed  lines  in  Fig.  11.  Now,  the  input  to  the  system 
is  the  controlled  relative  position  of  the  proof-mass  dc,  while  d 
represents  the  actual  relative  position.  The  accelerometer  and 


proximeter  feedback  gains  are  K  and  K  ,  while  A,  and  A„  are  control 

a  p  12 

gains.  For  K  =  1  the  equations  of  motion  appear  as 


Xj  =  (-KXj+  Ajind  +  Agmd  -  Ajmd^J/fM  +  m)  (81) 


d  =  -A.d  -  A0d  +  A,d 
X  »  1  c 


and  the  resulting  control  force  is  given  by 


•  »  ft 

u  =  -(d  +  Xj)m 


(82) 


(83) 


The  state  equation  for  the  expanded  system,  using  the  state  vector 
j 

s=[x^x^dd]  ,  becomes 


where 


-K/(M+m)  0  A1m/(M+m)  A2m/(M+m) 

0  0  0  1 


-A  m/(M+m) 


To  determine  the  gains  A^  and  Ag  in  the  position  loop,  consider  a 
transfer  function  d/d 


d/d  =  A,/(s  +  A„s  +  A.) 
c  1  v  2 


where 


A,  =  w  and  A  =  2f  u 

1  p  2  p  p 


Select  fp=  0.707  to  minimize  the  settling  time  of  the  actual  relative 
position  in  the  position  loop.  The  percentage  overshoot  of  the  system 
will  be  4.325  %.  To  take  the  overshoot  into  consideration  the  maximum 
rattlespace  should  be  reduced  to 


d  =  1.217  x  10  m 
max 


Assume  that  the  external  excitation  to  the  system  will  not  exceed  the 
control  force  level,  i.e.. 


co^lx,  |  <  w^d  =  u  /m 

1  1  'max  p  max  max 


-  -■>  v»Vw  ^jnyr^jfT-Vv  jr*jr*rj»  '>  ">  ■>'*>'”>  v1 ' 


’w.w.v.v.’^.- 


r.-TT.-v. 


where  u 


is  selected  from  Eq.  (89) 


u  /(md  )  or  u  =24.3  rad/s 
max  v  max7  p 


(90) 


Under  these  conditions  the  system  will  not  violate  the  given 
constraints.  Eq.  (74). 


The  input  d  is  determined  using  the  identification  method  based  on 
c 

fitting  optimal  LP/MST  trajectories,  i.e.,  d  (k)  and  s  (k)  for  k  = 

1.2 . N.  The  task  is  to  find  gain  d  such  that 


d*(k)  =*  d  s*(k)  =  dS(k) 


(91) 


where 


\  \  dd  dd  1 


(92) 


Note  that  d,  should  be  set  equal  to  zero  to  obtain  meaningful  curve 
d 

fitting  results.  Using  a  curve  fitting  technique,  suboptimal  feedback 
gains  are  obtained  as 


d  =  [-0.4439  -0.3878  0  -0.0349] 


(93) 


and  the  resulting  closed  loop  eigenvalues  are 


-31.40  +  26 . 89 i  and  -2.092  +  2.595i 


(94) 
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